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Abstract Ultraluminous X-ray sources (ULXs) are 
amongst the most intriguing of X-ray source classes. 
Their extreme luminosities - greater than 10 39 erg s _1 
in the 0.3 - 10 keV band alone - suggest either the pres- 
ence of black holes larger than those regularly encoun- 
tered in our own Galaxy (the Galactic centre excepted) , 
or sources apparently radiating well above the Edding- 
ton limit. We review the insights afforded us by studies 
of their X-ray emission, focussing on what this reveals 
about the underlying compact object. In particular, we 
discuss recent deep observations of ULXs by the XMM- 
Newton observatory, and how the unprecedented data 
quality provided by this mission is starting to discrim- 
inate between the different physical models for these 
extraordinary X-ray emitters. 
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■i-i . 1 Introduction 



The first imaging observations of galaxies beyond our 
local group were conducted by the Einstein observa- 
tory in the period 1978 - 1981, and revealed a some- 
what unexpected result: many galaxies hosted one 
or more extra-nuclear sources with X-ray luminosities 
well in excess of those typically obs erved in our ow n 
galaxy and its nearest neighbours (jFabbiand Il989l ). 
At the time it was rather presciently acknowledged 
that if these were individual X-ray sources, their X- 
ray luminosity was difficult to explain without invok- 
ing massive black holes, or super-Eddington emission 
(jFabbiano fc Trinchierj|l987t l. 
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Throughout the 1990s the ROSAT mission ob- 
served hundreds of near by galaxies, detecting many 
more of these sources (Colbert fe Mushotzkvl 1999: 



Roberts fc Warwidj200otlcolbert fc Ptalj200ilLiu fc Bregmarfr 
20051) . These observations revealed little of the actual 
nature of these objects, except that a fraction can be as- 
sociated with recent supernovae. for ex ample SN 1986J 
in NGC 891 (jBregman fc Pildisl 1 1 99^ ) . However, sub- 



sequent analysis of samples of these objects has given 
us some insights into their demographics, for example 
with estimates of between only 1 in 8 and 1 in 4 ma- 
jor galaxies hosting one of these X-ray bright objects 
with an observed luminosity in excess of 10 39 erg s _1 
jPtak fe Colbertll2004lLiu. Bregman fc Irwinll2006l ). It 
is this extreme luminosity threshold, in combi- 
nation with an extra-nuclear location, that we 
use to define this class of so-called "ultralumi- 
nous X-ray sources" (ULXs). 

A clearer insight into the nature of the majority of 
ULXs was provided by ASCA observations in the years 
around 2000. Most notably, the novel wide bandpass 
CCD spectroscopy afforded by ASCA allowed the spec- 
tra of many ULXs to be measured over the 0.5 — 10 
keV range for the first time. This revealed some of 
them to be well-fitted by the multi-colour disc black- 
body model, that describes the optically-thick thermal 
X-ray emission of an ac cretion disc around a black hole 
jMakishima et alj |2000) . Additionally, multi-epoch ob- 
servations of some ULXs saw them apparently tran- 
siting between spectral states described by cither the 
multi-colour disc blackbody model or a power-law, sim- 
ilar to the transition between low- and hi gh-states seen 



in G alactic black hole X-ray binaries (jKubota et al 
20011) . This gave the first strong corroborating evidence 
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- on top of their extraordinary luminosities - that most 
ULXs are accreting black holes. However, these results 
were not without proble ms, notably that t he dis c tem- 
peratures measured by Makishima et al. ( 2000h were 
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far too hot to be reconciled with the masses suggested 
by the high luminosities of the sources (see below). It 
was suggested that this could be explained by rapidly 
spinning (Kerr metric) black holes, in which the inner 
edge of the accretion discs are far closer to the black 
hole - and therefore hotter - than in the non-rotating 
(Schwarzschild) case. 

Our understanding of the nature of ULXs has in- 
creased immensely over the last ~ half decade due to 
the excellent capabilities of the Chandra and XMM- 
Newton observatories, and follow-up studies across the 
range of the electromagnetic spectrum. In this paper 
we mainly concentrate on contributions to the under- 
standing of ULXs garnered from spectroscopic and tim- 
ing studies using the European Photon Imaging Cam- 
era (EPIC) on XMM-Newton. However, first we will 
summarise some of the main arguments relating to the 
nature of ULXs. 

1.1 A new class of black holes? 

The Eddington limit for the maximum radiative lumi- 
nosity possible from the spherical accretion of matteiQ 
on to a black hole can be expressed as 



L 



Edd 



1.3 x lO d8 (M/M ) erg s" 



(1) 



where M is the mass of the acc reting object in so- 

Hence, for an ob- 



lar masses (Makishima et al 



ject obeying the Eddington limit, at a luminosity of 
10 39 erg s _1 its mass must be >7.7M©. For an ac- 
cretion rate of ~ 10 per cent of that required to reach 
the Eddington limit - a fairly typical accretion rate for 
a high-state black hole - this means that a ~ 77M© 
black hole is required for the source to be emitting at 
10 39 erg s _1 . So, if ULXs obey the Eddington 
limit, they must contain massive black holes. 
But how massive? 

Dynamical friction arguments imply that these 
sources cannot be misplaced super-massive black holes, 
sitting outside the nuclei of the host galaxies, as such 
massive objects sh ould sink to the centre of the galaxie s 
in a Hubble time (jTremaine. Ostriker fc Spitzerlll975l ). 
However, the Eddington limit argument also rules out 
the stellar remnant black holes that we know of in our 
own ga laxy, with masses in the range 3M@ < Mbh < 
18M (|McChntock fc Remillardll2006t ). for all but the 



mildest of ULXs (L x < 2.3 x 10 39 erg s" 1 at the Ed- 



dington limit). Indeed, iFrver fc Kalogeral (|200ll ) cal- 
culate that the vast majority of black holes formed 



1 This equation is strictly correct only for ionised hydrogen; the 
accretion of helium and/or heavier elements will raise this limit. 



from the evolution of a single massive star will have 
mass < 20-Mq, clearly inadequate to power the brighter 
end of the ULX population (if obeying the Edding- 
ton limit) iH These limits led to the suggestion that 
ULXs may be the first observational evidence for a 
new, ~ 10 2 — 10 5 Mq inte rmediate-mass class of accret - 
ing blac k holes flMBHsl IColbert fc Mushotzkvlll999l ) 
(see also lMiller fc Colbertl l|2004h for more on IMBHs). 

The strongest supporting evidence in favour of 
IMBHs in ULXs comes from the high signal-to-noise 
broad-band X-ray spectr oscopy enab l ed by XMM- 



Newtor^. In particular, Miller et al. ( 2003 ) showed 



that the spectra of two ULXs in NGC 1313 could be well 
fitted by the same absorbed multi-colour disc black- 
body plus power-law continuum model that is com- 
monly used as the empirical model to fit Galactic black 
hole binaries, with the key difference being a lower disc 
temperature in ULXs than Galactic black holes (0.1 
- 0.3 keV versus ~ 1 keV, respectively). This is cru- 
cial because, for a fixed accretion rate, the temper- 
ature of the inner edge of a standard accretion disc 
scales with the black hole mass as T oc M~ a - 25 , i.e. 
a cooler disc implies a bigger black hole. In fact, 
the black hole masses for the ULXs in NGC 1313 
were estimated to be of the order ~ 1000M Q . Many 
ULX s pectra were quickly shown to a gree with this 
result dMiller. Fabian fc Miller] l2004al: ICropper et al 



| 2004 iDewangan et alJl20q4 IRoberts et al.ll2005l ). with 
iMiller. Fabian, fc Miller] (|2004bl ) demonstrating that 
such sources lie in a different region of disc luminos- 
ity - disc temperature space than Galactic black holes, 
emphasizing their potentially different natures. 

Other factors also argue for the presence of IMBHs 
in at least some ULXs. For example, X-ray timing 
characteristics such as the detection of quasi-periodic 
oscillations (QPOs) in the X-ray fluctuation Power 
Spectral Densities (PSDs) of some ULXs argue that 
their emission is isotropic, supportive of IMBHs as- 
suming that the Eddington limit is not exceeded (see 
Section [3] for more details). Similarly, simple pho- 
ton counting arguments for high-excitation optical 
line emis sion regions near ULX s make the same ar - 
gument (jPakull fc Mirionil 120021 iKaaret et aD |2004 ). 
The source for which most evidence stacks up is M82 
X-l, which through a combination of its extreme lu- 
minosity (ix.pcak ~ 10 erg s _1 ), co-location with 



^Frver fc Kalogeral i200lh do note that it is possible that very 
massive, low metallicity stars leave a sufficiently massive rem- 
nant core after the wind mass-loss phase ( > 42Mq ) to collapse 
directly to a massive black hole. However such objects would be 
comparatively rare. 

3 Although similar results have been obt ained with ASCA 
dColbert fc Mushotzkvlll99Sh and Chandra jKaaret et alj|2003l : 
iRoberts fc Colber<Jl2003l: iRoberts et al.ll2004T ). 



X-ray observations of ultraluminous X-ray sources 



3 



the young, dense stellar cluster MGG 11, and QPO 
detections is the best known candidate for an IMBH 



(Kaaret et al. 2001: 


Strohmaver & Mushotzkv 


Portegies Zwart et al. 


2004; Mucciarelli et all 



12006|). 

However, it is possible this source is an atypical ULX; 
it may be the nucleu s of an accreted dwarf galaxy 
( King fc Dehnenll2005l ). 



1.2 The problem(s) with IMBHs 

Unfortunately, ULXs as a population are not triv- 
ially explained by the presence of IMBHs. There are 
in fact many arguments as to why the majority of 
ULXs cannot be IMBHs, or at least IMBHs of the 
size inferred from cool accretion discs (~ 1000M Q ). 
Two arguments stand out as the principle reasons 
the ULX population is not dominated by (large) ac- 
creting IMBHs. Firstly, the luminosity function of 
X-ray sources in galaxies (XLF) has an unbroken 
power-law form for 5 decades up to a luminosity of 
~ 2 x 10 40 erg s" 1 (|Grimm. Gilfanov fc Sunvaevll2003 : 
Swartz et al.l 120041 ). This break occurs at ~ 10 per 
cent of the Eddington luminosity for the ~ lOOOM© 
black holes inferred from ULX spectroscopy. This is ex- 
tremely troublesome for a ULX population dominated 
by these large IMBHs, as they would not only have to 
contrive to take over the XLF smoothly from Galactic 
black holes at ~ 10 39 erg s _1 , but then cease accret- 
ing at 10 per cent of Eddington. No other accreting 
source class behaves in this manner. This instead ar- 
gues that ULXs are dominated by black holes of mass 
up to ~ IOOA/q (or less if the Eddington limit can be 
exceeded). 

The second strong argument against IMBHs comes 
from the association of ULXs with star formation. 
Early observations with Chandra revealed that star- 
burst galaxies have populations of multiple ULXs 
( Fabbiano. Zezas fc Murravl200lllLira et al.l2002tlRoberts 
2002). an unusual result given the expectation of less 
than one in four galaxies on average possessing even 
one ULX. The obvious conclusion from this is that the 
ULXs are intrinsically linked to the ongoing star for- 
mation occuring in those galaxies. However, the direct 
co-location of ULXs with the star formation, most no- 



1.3 The most extreme stellar-mass black holes? 

If ULXs in starburst galaxies are not IMBHs, then 
what are they? The obvious solution is to turn to 
a class of objects we would expect to find there any- 
way: high-mass X-ray binaries. The problem then be- 
comes one of making such objects appear as ULXs. 
Assuming Galactic black hole masses for these objects 
(i.e. M < 2OM0) one then needs to either make 
such objects actually break the Eddington limit, or 
to make them apparently exceed it. Lower-luminosity 
and/or Eddington-limited objects could appear as 
bright ULXs due to beaming, either through relativistic 
boosting of their X-ray emission a long our line-of-sight 
(jKording. Falcke fc Markofj|2002j) or through collima- 
tion of their radiation - probably by a geometrically- 
thick accretion disc - such that it on ly escapes into a 
fraction of the sky ( King et al.ll200lM 4 l Alternatively, 



models have been suggested whereby actual super- 
Eddington luminosities are achieved and maintained, 
at factors < 10 above the Eddington limit ( Bcgclman 



2002tlEbisawa et alJliooa iHeinzeller fc Dusch]|200 



Regardless of the processes involved, a very basic 
requirement of most models is that sufficient fuel is 
available for the s uper-Eddington mass transfer rates 
neede d in ULXs. iRappaport. Podsiadlowski fc Pfahll 



(2005) show that this is indeed the case for high-mass X- 
ray binary systems containing a stellar-mass black hole 
and a massive donor star, that can in fact sustain super- 
Eddington mass transfer over a very large fraction of 
their lifetimes. Other authors suggest that the possi- 



objects fuel ULXs 


Begelman, King & Pringle 


2006; 


Poutanen et alJl2007 


. Furthermore, where optical stel- 



lar counterparts to ULXs have been identified, primar- 
ily by HST, they tend to be blue and of an appropriate 
magn it ude for the young, ma ss ive stars required to fuel 
WM M dRoberts et al.l |2001[ iLiu. Bregman fc Seitzerl 



2004lKuntz et alj2 005) . This provides compelling sup- 
port to the argument that ULXs arc high-mass X-ray 
binaries. 

Finally, it is reassuring to know that the Edding- 
t on limit is broken in prac t ise in Galactic black holes 
McClintock fc Remillardl ( 2006 ) give several exam- 



tably seen in the Cartwheel galaxy (|Gao et al.l 12003). 
implies that they must be (relatively) short-lived, which 
requires successive generations of ULXs to be fo rmed 
over the duration of the star formation event. King 



( 20041 ) pointed out that if these ULXs were all large 
IMBHs, then an infeasibly large proportion of the avail- 
able star forming mass would end up in the form of 
IMBHs. Hence the majority of ULXs in star forming 
regions cannot be powered by IMBHs. 



pies, most notably that of GRS 1915+105. This source 
is sufficiently luminous (at ~ 10 39 erg s _1 ) to ap- 
pear as a ULX if viewed from outside our galaxy, 



4 QPO detections and high-excitation optical line measurements 
argue against all but the mildest forms of beaming in some ULXs. 
Additionally, the lack of detectable radio emission and/or rapid, 
high amplitude X-ray variabili ty in most ULXs arg ues against rel- 
ativistic beaming - though see iKrauss et ail l|2005l 'l for a possible 
counter-example. 
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and (for its known black hole mass) has consistently 
displayed peak luminosities in excess of the Edding- 
ton limit over the ~ 15 years of its outburst to date 
( Done. Wardzihski fc Gierlihskill2004l ). As the Edding- 



ton limit is exceeded in known sources, there is no rea- 
son this cannot also be occuring in ULXs. 

We are therefore left with the situation where most 
ULXs could be explained by stellar-mass black holes 
that are either super-Eddington, or subject to some sort 
of beaming. However, it is still difficult to reconcile the 
most extreme ULXs - those above 10 40 erg s _1 - with 
simple stellar-mass black hole systems. Larger black 
holes would still provide an obvious solution. But are 
they really the ~ 1000M Q black holes inferred from 
ULX spectroscopy? 



2 X-ray spectroscopy 

2.1 The IMBH model, and other solutions 

As we have already discussed, the strongest support 
for ULXs containing IMBHs comes from X-ray spec- 
troscopy, and in particular the good fits obtained to 
ULX spectra using the same multi-colour disc plus 
power-law continuum model used for Galactic black 
holes. To date more than 10 ULXs with decent XMM- 
Newton spectra have been shown to have their spectral 
fits improved substantially (compared to, say, a simple 
absorbed power-law model) by the use of this model, 
in many cases producing a statistically acceptable fit 
to the data. An ubiquitous feature of these fits is a 
cool disc which, as we have already seen, when inter- 
preted at face value implies the presence of an IMBH 
with mass of the order ~ 1000M©. 

But should we take this mass at face value? There is 
good reason not to. The mass estimates deriving from 
the multi-colour disc blackbody model assume that this 
is the dominant emission component in the X-ray spec- 
trum (essentially, that the black hole is in the "high" 
state). Unfortunately, in ULXs fit by this "IMBH 
model" (which here specifically refers to an empirical 
model composed of a cool disc plus harder power-law) 
it is very evidently not. We demonstrate this graphi- 



cally for the IMBH candidate (L x > 10 40 erg s 



- 1 ) NGC 



1313 X-l in Fig[TJ where it is obvious that the power-law 
dominates the flux within the XMM-Newton bandpass. 
In fact, the disc typically emits no more than 20 per 
cent of the 0.3 - 10 keV flux of U LXs in this model 
( Stobbart. Roberts fc Wilms! l2006h . Worse still, the 
power-law slopes measured by this model are somewhat 
on the low side for the classic high state (r ~ 1.6 — 2.5 
for IMBH models, compared to T ~ 2.1 — 4.8 in Galac- 
tic high state sources; cf. McClintock fc Remillardl 




1 2 
Energy (keV) 



Fig. I.— XMM-Newton EPIC-pn spectrum of NGC 
1313 X-l, obtained on 2000 October 17, and re- 
processed with SAS version 6.5.0. The data points are 
shown in black, with the best-fitting IMBH model in 
green. The contributions to this model of the multi- 
colour disc black body component (with kT- ln ~ 0.2 
keV) and the power-law continuum (F ~ 1.7) are shown 
by blue (dashed) and red lines respectively. The X- 
ray emission detected by XMM-Newton is clearly dom- 
inated by the power-law component, and not the disc. 



(2006)). This makes black hole masses obtained by this 
method (at the very least) questionable. 

Furthermore, the IMBH model is not the only model 
that fits ULX spectra. Several sources have been iden- 
tified in which a variant of this model, where the disc 
component fits to the h ard end of the spectrum, pro- 
vides a far superior fit (IStobbart. Roberts fc Warwick 



20041: iFoschini et"aI1l2004 Feng fc Kaaretl 120051) . As 



discussed by Roberts et al. ( 20051 ) this variant of the 
empirical black hole spectrum model does not provide 
a physical model for the X-ray emission - for exam- 
ple there cannot be sufficient photons present in the 
vicinity of the black hole to produce the dominant soft 
power-law through Compton up-scattering. However, 
the crucial point is that there is distinct curvature 
(which can also be described as a spectral "break") 
present above 2 keV in these data, for which there must 
be a physical explanation. 

2.2 Re-evaluating ULX spectra from XMM-Newton 

Given the question marks about the IMBH model, and 
this second, "inverted" model that fit to some ULXs, 
we set out to examine the best available XMM-Newton 
ULX d atasets (from the beginning o f 200 5). Specifi- 
cally, in IStobbart, Roberts fc Wilmsl (|2006h we set out 
to ask the questions 
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(a) Power— law fit to pn data, ULX in M83 




1 2 
Energy (keV) 

Power-law fit to pn data, M33 X-8 
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Fig. 2 — XMM-Newton EPIC-pn data fo r two ULXs 
in the IStobbart. Roberts fc Wilms! (|2006!) sample. In 
both cases we show the data (in red) and best-fitting 
power-law continuum model (in blue) in the top panel, 
with the A\ residuals for the fit shown in the bottom. 
The ULX in M83 was adequately fit by a power-law 
continuum, whereas M33 X-8 was not. 



• How easy is it to distinguish the IMBH and in- 
verted models given the available quality of data? 

• With what frequency do these models work for 
the available data? 

• Can we say anything about the physics of the in- 
verted model? 

To do this we selected data sets with at least a few 
thousand counts (EPIC-pn and MOS combined) per 
ULX. We show examples of the low and high end of the 
data quality for ULXs in our sample in Fig [2] The data 
was of sufficient quality to statistically rule out simple 
multi-colour disc blackbody fits to all the sources, and 
power-law continua in 8/13 cases. 

As simple models were inadequate for all but 
the poorest data, we next attempted empirical two- 



component models, beginning with the IMBH model. 
This provided acceptable fits to 8/13 sources, which all 
had the classic ~ 0.1 — 0.3 keV cool disc signature. Un- 
fortunately, they all also displayed the problems inher- 
ent in mass measurements from this model, i.e. domi- 
nant, hard power-law continua. We then attempted fits 
using the inverted model. This also provided good fits 
to 8/13 data sets, parameterised by T ~ 2.5—4.3 power- 
law photon indices, kT ln ~ 0.9 — 2.7 keV inner-disc tem- 
peratures, and a very roughly 50/50 split between the 
flux contribution of the two components in the XMM- 
Newton band. Interesting, the six lowest quality data 
sets provided acceptable fits to both models, demon- 
strating that either very high quality X-ray data, or 
a secondary diagnostic, is required to distinguish the 
models. 

This key diagnostic is found at energies above 2 keV, 
where disc-domination leads to distinct curvature in the 
spectrum, whereas a dominant power-law has an unbro- 
ken spectrum. We tested for this characteristic signa- 
ture by comparing power-law fits to broken power-law 
fits on the > 2 keV data for each ULX. In total we found 
8 ULXs showing evidence for breaks, at significance lev- 
els between 3 — 10ct. Three of these were expected, as 
they were from sources clearly better fit by the inverted 
model. However, five sources that were either ambigu- 
ous or well-fitted by the IMBH model also showed a sig- 
nificant break. Of the remaining (unbroken) source fits, 
in three cases this may be attributable to very poor data 
quality above 2 keV. This simple test therefore demon- 
strates that spectral breaks are present in the majority 
of ULXs (notably, across the whole range of luminos- 
ity) where the data quality above 2 keV is sufficient to 
detect themU 

In order to investigate the physics of this break 
further, we attempted spectral fits using a physically 
self-consistent accretion disc plus Comptonised corona 
model, specifically using an absorbed DISKPN + EQ- 
PAIR model in xspec. Where appropriate, parame- 
ters were tied to assumed values (based on experience 
with Galactic systems), so that the model had only 
one more degree of freedom than the empirical (two- 
component) models, but sufficient scope for variation 
within the model parameters was allowed such that the 
outcomes were not pre-judged. This model gave supe- 
rior fits to the empirical models, with 11/13 ULXs pro- 
viding statistically-acceptable fits (with one more only 

A by-product of this spectral break/curvature is that the 
best empirical fits are provided by models composed of 
two thermal components, in particular a combination of a 
blackbody with a multi-colour disc blackbody model - see 
IStobbart, Roberts fc Wilms! (|2006h for more details. 
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marginally unacceptable). Two remarkable character- 
istics were common to most fits: firstly, the ULXs still 
showed apparently cool discs; but secondly, the coronae 
(in 9/12 cases) appeared optically-thick (with t > 8 
- in several cases much higher). It is this optical thick- 
ness in the corona that is responsible for the curvature 
in the 2-10 keV spectrum. Furthermore, this distin- 
guishes ULXs from galactic black holes, that do not 
typically show such coronae. The reason their coronae 
are empirically modelled by power-laws is that they are 
optically-thin. This therefore suggests ULXs are oper- 
ating in a different accretion mode to the classic states 
of Galactic black holes. 

2.3 Possible scenarios 

So what is happening in these sources? One physical 
scenario that could lead to the production of both a cool 
disc com ponent and an opt ically-thick corona is sug- 
gested bv lZhang et al. ( 2000h . based on an analogy with 
the solar corona. This "sandwich" model consists of a 
cool inner accretion disc, seeding an outer, warm (~ 1 
keV) and optically-thick (t ~ 10) accretion disc layer 
with ~ 0.2 keV photons. This could readily provide 
both components we derive from our spectra, though 
the geometry required to see the cool photons through 
the warm layer is problematic. A second scenario is 
based on observations of a strongly-Comptonised very 
high state ob served in the Galactic black hole XTE 
J1550-564 bv bone fc Kubotal (|2006l ). In this source 
they also detect a cooler-than-expected disc, alongside 
an optically thick corona (though neither phenomena 
are as extreme as detected in the ULXs), and suggest 
that this is due to energetic coupling of the inner-disc 
with the corona. In such a system, the energy released 
by the extreme accretion rate is sufficient to launch an 
optically-thick corona, which obscures (and, through 
extracting the launch energy, cools) the central regions 
of the disc. One therefore predominantly sees the cooler 
outer regions of the accretion disc, in addition to the 
optically-thick corona. As this state occurs at a very 
high accretion rate in XTE J1550-564, this suggests 
that ULXs operate at similar - or higher - accretion 
rates. This must be at around the Eddington limit, 
suggestive of black hole masses up to ~ 100 Af Q . 

Further work has now revealed this > 2 keV 
break in other ULXs - for exampl e in M82 X-l 
( Okaiima. Ebisawa fc Kawaguchi 2006f) and Ho IX X - 
1 faka M81 X-9) (|Dewangan. Griffiths fc RaollioOrj ). 
Other ideas have been also been postulated for its ori- 
gin, and the full observed X-ray spectrum for ULXs. 
Several authors have discussed ULXs in the context o f 
slim disc models (jWatarai. Mizuno fc Mineshigdbo'oih . 



In particular, they note that as the Eddington limit is 
approached, the structure of the accretion disc should 
change. This would manifest itself as a change in the 
model disc profile, T(r) oc r~ p , where standard discs 
have p = 0.75, and slim discs p = 0.5. Recent work 
where a variable disc profile model is fit to ULX spectra 
does i ndeed show values of p ~ 0.6, suggestive of slim 
discs |Vierdavanti et al-lEoO^ iMizuno et alJ l2007h. A 



second idea, put forward by Goncalves fc Soria ( 2006h . 
draws from observations of AGNs with outflows. In 
their model ULXs have an intrinsic very high state 
spectrum (steep power-law form, i.e. T > 2.5), that is 
modified by absorption from material in an ionised fast 
outflow. This effectively takes a "bite" out of the spec- 
trum in the ~ 1 — 4 keV range, resulting in the apparent 
soft excess and > 2 keV spectral break features. We 
note that the post-break power-law slopes measured by 
Stobbart. Roberts fc Wilms! (|2006l ) are indeed consis- 
tent with a very high state spectrum. Finally, work by 
Poutanen et al.1 (|2007N has related ULX spectra to the 
expected spectra from super-critically accreting black 
holes viewed with the accretion disc close to face-on. 
In such sources - SS 433 may be an example, viewed 
edge-on - the break originates from a direct view of the 
spectrum of the inner-regions of a hot accretion disc. 

All these models have a common theme - accretion 
at around or above the Eddington rate. This strongly 
suggests black holes of <100M Q , rather than the ~ 
1000M Q IMBHs formerly proposed. 

2.4 Spectral variability 

Snapshot observations, though useful, can only tell 
us so much about accreting sources. As 10 years 
of RXTE observations have shown, for example with 
Galactic black holes, further progress can be made by 
consid ering how these observed spe ctra change with 
time |McClintock fc Remillardl l2006h . As ULXs are 
in general too X-ray faint for monitoring missions 
like RXTE (with the excep tion of M82 X-l - see 
Kaaret. Simet fc Lang! (|2006l )1. the best one can do is 
to use monitoring campaigns of snapshot observations 
separated by days - months. 

Unfortunately, to date very few campaigns have 
been pursued. Where they have been undertaken, 
sources are observed to behave in two distinct ways. 
Some ULXs behave as would be expected from clas- 
sic Galactic black hole behaviour, i.e. they get spec- 
trally softer as their flux increases. In contrast, 
other U LXs behave in the opposite fashion - they get 
harder ( Fabbiano et al. _ 2003t Dewangan et al. _ 2004 



Jenkins et al.l2004tlFeng fc KaaretJ2006al lbl: ISoria et al. 
200/1 We have undertaken such a campaign for the 
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Fig. 3. — Chandra ACIS-S hardness ratio - count rate 
plots for NGC 5204 X-l. The bands used are S = 0.3-1 
keV; M = 1 - 2 keV; H = 2-8 keV. Blue data points 
are from an initial 50-ks observation, with data from 
ten 5-ks follow-ups shown in red. The ULX spectrum 
clearly hardens as its flux increases. 



ULX NGC 5204 X-l using Chandra data, with ob- 
servations separated by days - weeks, and find that 



it behaves in the latter mode (jRoberts et al.l [2006) 



see Fig [3] In particular, we show that this behaviour 
can be modelled as changes in the temperature of the 
optically-thick component of a cool accretion disc plus 
thick corona model, with the corona heating up as the 
luminosity of the ULX increases. This behaviour was 
also seen in the str ongly-Comptonised very high stat e 
of XTE J1550-564 (jKubota fe Done! 120041 : fsoriall2007h . 
confirming the viability of the optically-thick corona 
model for this ULX. 



3 X-ray timing 

Perhaps the most fundamental timing signal from a 
ULX is an X-ray periodicity, as it can be used as a 
first step towards establishing the orbital characteris- 
tics of the underlying binary system. However, such 
measurements are not very common, with the best 
detections coming in a couple of cases where eclipses 



have been found ( Bauer et al.l l2001t iDavid et all 12005 



Fabbiano et al. 20061 ). Though other claims have been 
made, they suffer from a lack of data - most are based 
on <C 10 cycles, and require confirmation through fur- 
ther observations. In the absence of periodicities, the 
most useful timing diagnostic is Power Spectral Density 
(PSD) measurements for ULXs. 



Fig. 4. — A simplified linear relationship between black 
hole mass and PSD break frequency for Cygnus X-l and 
a number of AGN. This plot is a reasonable approxi- 



mation for the Tb re a 



BH 



'Edd 



relationship found 

" if JJll E/UU * 

( 20061 ). assuming sources at the Ed- 
dington limit (where Tbroak is the break timescale, Mbh 
is the black hole mass, and rriEdd is the accretion rate in 
Eddington units). Note that sub-Eddington accretion 
rates move the break to the left on this plot (and vice 
versa). We shade the region for which XMM-Newton 
observations are sensitive, showing that this coincides 
with IMBH-like masses. Courtesy S. Vaughan. 

3.1 Power Spectral Densities for ULXs 

The shape of the fluctuation PSD for ULXs is poten- 
tially a very powerful tool for ULXs. In particular, 
the characteristic frequency of breaks in the PSD slope 
can be used to infer masses based on a direct scaling 
of properties betwe e n Gal actic black holes and AGN - 



McHardv et al.l (|2006l ). In Fig H] we demonstrate 



see 

this concept using a simple linear scaling of black hole 
mass to break frequency, and show that the 10 - 1000 
s timescales over which XMM-Newton observations are 
most sensitive to measuring breaks matches well with 
that in which one might find IMBHs. 

Unfortunately, ULX s do not generally show much 
short-term variability (jSwartz et al J 120041 ). with very 
few examples displayi ng sufficient variabilit y power to 
establish a PSD from (jFeng fc Kaaretj|2005l). However, 
some measurements have been made! ICropper et al.l 
( 20041 ) detected a putative break at 28 mHz in the 
PSD of NGC 4559 X7, which they suggest supports 
the case that it harbours a ~ 1000M Q black hole, al- 
though tjie_rjres^nc^ofJhis spectral break is now dis- 



puted ([Barnard et al.l 120071 ) . Another detection of a 



break frequency was made by ISoria et al.l (|2004h . who 
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found a break at 2.5 mHz in the PSD of NGC 5408 X- 
1. When the shape of the PSD either side of the break 
was considered, the a uthors derived a mass of ~ lOOM^ 
for thi s ULX. Finally, iDewangan. Titarchuk fc Griffiths 
|2006l ) find a break at - 34 mHz in the PSD of M82 
X-l, from which they infer a mass in the 25 — 520M Q 
range for the underlying black hole. 

In 2004 we were awarded a 100-ks XMM-Newton ob- 
servation in order to derive a PSD for the nearb y, lumi- 
nous ULX Holmberg II X-l (jGoad et al.ll2006l ). How- 
ever, we found that it displayed no strong variability 
during the observation. In fact, through a PSD anal- 
ysis we were able to demonstrate that the fractional 
variability amounted to less than a few per cent rms 
over time scales of minutes to hours. This variability 
power is less than that observed in classic high and very 
high state (red noise) PSDs. The PSD could instead be 
consistent with the low/hard state, but the energy spec- 
trum for Ho II X-l is closer to a very high state spec- 
trum (in fact, it displays a cool disc plus optically-thick 
corona spectrum) . One solution is that the source could 
be in a state similar to the 'x'-class of GRS1915+105, in 
which the source is in a very high state, but its PSD is 
band-limited (i.e. all its variability power is limited to 
a narrow frequency band) . As we do not see this band- 
limited variability, it is likely to be at higher frequencies 
than we are sensitive to in our PSD. This means the 
black hole must be small. In fact, we calculate a limit 
of < 100M Q for Ho II X-l from the lack of variability. 

3.2 QPOs 

A second feature of PSDs that has diagnostic potential 
for ULXs are quasi-periodic oscillations (QPOs). In 
fact, the detection of 3:2 ratio twin-peak high-frequency 
QPOs at ~ 1 Hz wo uld be strong evidence for the pres- 
ence of an IMBH |Abramowicz et all liooi ) (though 
scaling from Galactic systems implies such a measure- 
ment is unlikely even in the medium term). Low fre- 
quency QPOs, on the other hand, have now been de- 
tected in a handful of ULXs ( Strohmavcr & Mushotzkv 



20031 IDewangan. Griffiths fc Raoll2006l : lMucciarelli et al.| 
200(jt Istrohmaver et al.ll2007t ). However, this type of 
QPO docs not provide a clear, unambiguous mass 
estimate. For example, the detection of a QPO in 
NGC 5408 X-l, even with the additional diagnostics 
of a break frequency and a second (4:3 ratio) possi- 
ble QPO, could imply mass estimates anywhere in the 
range 100 — lOOOM^ depende nt upon the assumptions 



made ( Strohmaver et al 



20071 ). Despite this, QPOs are 



telling us one thing - as they are a coherent signal, their 
detection rules out all but the mildest forms of beaming 
in ULXs where they are present. 



4 Concluding remarks 

New observational evidence is now pointing away from 
the interpretation of ULXs as the ~ lOOOA/© black 
holes inferred from simple cool disc plus power-law 
spectral models. Putting to one side the inherent prob- 
lems with the dominant, hard power-law component 
derived from this model, the crucial evidence in this 
matter is the spectral break abov e 2 keV detected by 
Stobbart. Roberts fc Wilms! (|200d ). This is completely 
unexpected and inexplicable in the context of the simple 
~ 1000M Q black hole interpretation for ULXs. Indeed, 
it suggests physical characteristics that have more in 
common with Galactic sources accreting at around the 
Eddington limit. This implies much smaller black hole 
masses - of the order ~ 100M Q or less - for ULXs. 

Are ULXs then stellar-mass black holes that are ra- 
diating at super-Eddington rates, i.e. factors of < 10 
above the Eddington limit for most ULXs, assuming 
masses of < 20M Q ? Though this is feasible, it is per- 
haps not necessary. Almost all ULXs could trivially be 
explained (at least in luminosity terms) by accreting 
black holes with masses of a few tens of M Q , consis- 
tent with their spectra suggesting they are accreting at 
around (including slightly above) the Eddington limit. 
Interestingly, results based on optical observations of 
ULX counterparts are suggesting a similar conclusion. 
For example, a possible radial velocity variation in a 
He II line de tected from the optical s pectr um of NGC 
1313 X-2 bv IPakull. Grise fc Motchl (fcood) suggests a 
small black hole mass. Furthermore, irradiation mod- 
els of donor stars in ULXs, when combined with op- 
tical colou rs, suggest black hole m asses of the order 
<100M Q ijCopperwheat et alJl2007h . A final piece of 
the puzzle is that the creation of such black holes may 
be possible in the young stellar populations that we 
generally find ULXs co-habiting with, from the merging 
of a binary composed o f very massive early-type stars 
( Belczvnski et al. 2006) (or, alternatively, see footnote 
2). This could yield the black holes with masses of up 
to ~ 100M© that we are potentially finding in ULXs. 

Hence it now appears that we can tentatively con- 
clude that the vast majority of ULXs could harbour 
the slightly bigger cousins of the Galactic black holes 
we are familiar with, rather than their far larger dis- 
tant relatives. Should we call them IMBHs? The def- 
inition of this term is somewhat indistinct, with some 
authors quoting 20M Q as the lower limit for this class, 
and others starting at ~ 100Af Q . The few times larger 
than stellar-mass black holes that could be powering 
ULXs are precisely in this grey area. Perhaps a new 
artificial distinction between "small" ( < 100M Q ) and 
"large" IMBHs is required for clarity's sake. 
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However, the mass of an underlying black hole is 
yet to be conclusively determined for any individual 
ULX - this requires the measurement of a dynamical 
mass function, which is non-trivial for extra-galactic 
sources. Furthermore, a small and very rare sub-group 
of the most luminous ULXs - the "h yperluminous X- 
ray sources" with Lx > 10 41 erg s 



dGao et al.ll200i 



Wolter. Trinchieri fc Colpl boOfj iMiniutti et"aIT f2006) 



- defy easy explanation by anything other than large 
IMBHs, assuming they are indeed accreting black holes 
in the host galaxies, and not luminous supernovae or 
foreground/background objects. The issue of whether 
some ULXs could still constitute evidence for accretion 
onto large IMBHs is therefore far from finished with. 

TPR would like to thank his collaborators on the 
work presented here - most notably Ann-Marie Stob- 
bart and Mike Goad - and a select yet too-numerous- 
to-mention group of workers in this and related fields 
for many very useful conversations. He would also like 
to apologise for omitting to mention many other good 
pieces of work on ULXs due to simple space and sub- 
ject limitations. Many of the results quoted in this 
work are based on observations obtained with XMM- 
Newton, an ESA science mission with instruments and 
contributions directly funded by ESA Member States 
and NASA. Finally, thanks to Martin Ward for present- 
ing this work at the 5th Stromlo Symposium in TPR's 
absence. 
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